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1. Introduction 

The mixing of neutral B mesons provides a fertile area for tests of the Standard Model and 
sensitivity to new physics. In order to take advantage of experimental measurements, lattice com- 
putations are as usual needed to determine the hadronic matrix elements of the effective weak 
operators. The mixing is dominated by the short distance contributions from the operators [1] 

^1 = [by^Lq) {hfLq\ , G2 = [bLq) [bLq] , ^3 = (bLq] [bq) , 

^4 = (bLq) [bRq] , ^5 = {bLq] [bRq) , (1.1) 

where t&i through ^3 control the mixing in the Standard Model, while 1^4 and ^5 can appear in 
beyond-the-Standard Model (BSM) theories. Pairs of round or square parentheses in Eq. (1.1) 
indicate how the color indices are contracted, and R and L are the right and left projectors. Using 
Fierz transformations and parity, the mixing matrix element of any 4-quark operator with these 
quantum numbers can be written in terms of the matrix elements of these five "basis" operators; 
see, for example, Ref. [2] for details. 

In a lattice computation, it is useful to be able to fit the lattice data to a version of chiral pertur- 
bation theory that includes the effects of the discretization errors associated with the lattice action. 
In two recent lattice calculations of B mixing [3,4], staggered light quarks are combined with non- 
staggered heavy quarks using NRQCD [5] or the Fermilab action [6]. In such cases, the appropriate 
chiral theory is "rooted, heavy-meson staggered chiral perturbation theory" (rHMS^PT) [7]. Here, 
I describe a calculation of B mixing to one-loop order in rHMSZPT. Roughly speaking, the cal- 
culation is to leading order in the heavy-quark expansion, although the largest l/mg effects (the 
B-B* hyperfine splitting A* and the B^-Bd flavor splitting) are also included. This is a systematic 
approximation in the power counting introduced by Boyd and Grinstein [8] and discussed recently 
in Ref. [9] for the lattice calculation of heavy-light meson decay constants. 

The rHMSZPT calculation needs to take into account the form of the lattice operator used 
to approximate the continuum one. References [3, 4] construct the 4-quark operators as the local 
product of two local bilinears, each formed from a heavy antiquark field and a naive light-quark 
field. Both the use of naive fields and the local nature of the 4-quark operator influence the form of 
the coiTcctions at one-loop. 

2. Detailed Structure of the Operators 

Heavy-light bilinears and 4-quark operators are made by converting a one-component stag- 
gered fermion x{x) into a naive fermion ^{x) following Ref. [10], and then coupling it locally to 
the heavy-quark field Q{x): 

bilinear: Q{x)rW{x) , (2.1) 

4-quark operator: Q{x)r^{x) Q{x)r'^'{x) , (2.2) 

where F and F' are Dirac spin matrices. As we will see, bilinears of this form work exactly as 
desired in the lattice simulation, but the local product of the bilinears in the 4-quark operators 
introduces contributions from operators with wrong taste and spin (i.e., spin different from F(8)F'). 
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The naive light-quark action can be written as four copies of the staggered action: 

^>{x) = aix) X , ^{x) = r^" 7[' . (2.3) 

Here ;^ is a "copied" staggered field, with each Dirac component Xi separately having the staggered 
action. I will call the SU (4) symmetry that acts on index / "copy symmetry." Copy symmetry is an 
exact lattice symmetry, so copied and uncopied propagators are related by 

{Ux)piy)) = 5u'{xix)xiy)) ■ (2.4) 

where x is the normal (uncopied) staggered field. An interpolating field J^{x) for a heavy-light 
pseudoscalar meson is 

W = Q{x) 75 ^(x) = Q{x) 75 Q.[x) x{x) • (2.5) 

In the simulations, M'{x) is always summed over a time-slice (either explicitly, or implicitly by 
using translation invariance). To leading order in a, Q{x) varies smoothly between neighboring 
spatial sites (up to gauge transformation), but x does not, due to taste doubling. Staggered fields 
are, however, smooth in the spin-taste basis on the doubled lattice, so we need to sum the field 
within a hypercube to expose the structure of the operators. 

I focus first on the average of M'{x) over a spatial cube. Let x= {t,x) with x = 2y even, and 
let T] = (t]o> be a 4- vector with all components or 1. For t even (t = 2t), 




^^Qit,x)Y5Y.H'l)^(^^ + 'lo,2y + ri) + {-\)'^°a{T^)x{2r + T^o,2y + ri)] . (2.6) 

ID ^ 

Inserted gauge links for point-split quantities are implicit. For t odd, the result in Eq. (2.6) is the 
same except the last term changes sign. This is the usual oscillating state with opposite parity. 

For simplicity, I assume from now on that the oscillating state is removed by the fitting proce- 
dure, and that all components of x are even. Then 

^ (2(x)75 Z^{Tl)x{2y + 7l) ■ (2.7) 
16 ^ - 

To convert to a spin-taste basis, one can use [1 1] 

^r(3') = ^I^^""('7)^K2j + r7), (2.8) 

° ri 

where a is a spin index, a is a taste index, and / is a (trivially inserted) copy index. This implies 
that taste and copy indices are coupled in J^. We get (with spin indices implicit from now on) 

J^^'''\x)^^-Q{x)Y5qKy)8^' , (2.9) 

where the arrow signifies that oscillating states and i^(a) corrections are being dropped. Using 
Eq. (2.4), we see that the contraction of with Jif^ is automatically averaged over tastes: 

(jr(x) ^ 1 {Q(^)y,g"(y) q\y')y,Q{x')) . (2.10) 
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Note that the interpolating field constructed from the naive light quark gives us the desired pseu- 
doscalar spin in the spin-taste basis (aside from the dropped oscillating state). 

I now turn to the 4-quark operators constructed as in Eq. (2.2). The two bilinears are not 
separately summed over space, but they can be disentangled using 

^Y.^r{£l{^)K^\^)K) tr{a{r^')KQ.^ri')K) = 5,,,^. , (2.11) 

K 

where K runs over the 16 independent Hermitian gamma matrices. Following the same kinds of 
manipulations as in Eqs. (2.5) through (2.9), one then finds 

^n'"^ (x) ^ \ L(<2r„^^^, QKKqf) K^kKdi , (2.12) 

where r„ and are the Dirac spin matrices of the operators. Here the contributions with K ^ I 
clearly have the wrong (undesired) spin, as well as non-trivial coupling of taste {c,d) and copy 
{k,i) indices. Putting Eq. (2.9) together with Eq. (2.12) and using copy symmetry gives 

(^Vi"''^ ^) oc {D"'D''')KcaKde + (second equivalent contraction) , (2.13) 

where D"'^ is the Ught-quark propagator (in a given background) for taste a into taste c. When 
taste symmetry is exact, {D"'^D'^'^) oc deeded, so K = I, and only the coiTcct spins contribute. This 
shows that 4-quark operators constructed as in Eq. (2.2) will have the con^ect (desired) spin in the 
continuum limit, as well as at tree level in rHMS^PT, which respects taste symmetry. But, at one 
loop (and a / 0), taste-violations imply that {D"'^D^'^) is not necessarily proportional to SacSgd. 
Thus one-loop diagrams will have contributions from wrong-spin operators. 

3. Calculation of the ZPT Diagrams at One Loop 

We now are ready to perform the calculation of the one-loop chiral corrections to matrix el- 
ements of the 4-quark operators. Aside from the wave-function renormalization diagram, which 
does not involve the operators and can be taken unchanged from the literature [7, 12], there are 
three diagrams, shown in Fig. 1. 

Because copy and taste indices are coupled in Eq. (2.12), and because copy indices follow 
the quark contractions (Eq. (2.4)), we need to consider the quai^k flow of the meson diagrams 
in order to compute them. As an example. Fig. 2 shows a possible quark flow for each of the 
tadpole meson diagrams. These particular quark flows have "connected" pion propagators without 
hairpin vertices; other tadpole flows with either taste-singlet (physical) or taste-violating hairpin 
contributions are also possible. 

Labeling the taste of the pion in the loop by S (S = 1, . . . , 16), diagram Fig. 2(a) is propor- 
tional to 'Eaf'^fc^edKcaKde = [tr(K)]^, whcrc the factors of K come from Eq. (2.12). In this case, 
only the coiTcct spin contributes {K = I). On the other hand, diagram Fig. 2(b) is proportional to 
"Ead '^ec-KcciKde = tr(S K E K). If the propagator were independent of S, the sum on tastes would 
again give K = I. Due to taste violations, however, the pion propagator depends on whether E is P, 
A, T, V, or I. This implies that various K values are possible, giving wrong-spin contributions. 



4 



Staggered Chiral Perturbation Theory for Neutral B Mixing 



C. Bernard 




B B* B* B 



(C) 

Figure 1: Meson-level chiral diagrams for the mixing matrix elements. The gap between the two B (or two 
B*) fields represents the insertion of the 4-quark operator. Diagrams (a) and (b) are tadpole diagrams, which 
are distinguished by whether the contracting pions come from the same meson field (loosely speaking, the 
same bilinear) or from different meson fields. Diagram (c) is the sunset diagram. 




Figure 2: Examples of quark-flow tadpole diagrams. The filled, black circles represent the heavy and light 
quarks in the 4-quark operator. Diagram (a) contributes to the meson diagram Fig. 1(a), and diagram (b) 
contributes to diagram Fig. 1(b). Indices a-f label tastes. 



The fact that wrong spins enter means that other operators appear, and they in turn have dif- 
ferent chiral representatives than the original operator does. Fortunately, the basis of , . . . , ^5 is 
complete, and the chiral representatives of all these operators are given by Detmold and Lin [13]. 
Generalizing to operators with light-quark tastes c, d, we have 
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orcO(i], (3.1) 
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x,d 



x.d 



[orcf^J]. (3.3) 



where P and P* are heavy-light meson fields, a is the pion field, and x is the light flavor. The effect 
of copy indices is to enforce the contraction of the external quark of taste a with the quark in the 
operator of taste c, and similarly for e and d. 

We then write the matrix element for operator ^„ as 



{t.mB'D =^Ji + Wb + ^x^"^ + ^x'"' ) + jS; ( + ) + analytic terms 



(3.4) 



where #b is the B wave-function renormalization, ^ and are the right- and wrong-spin tadpole 
diagrams, and £2 and are the right- and wrong-spin sunset diagrams. In the special case of 
operator ^1, jSj = jSi by a heavy-quark spin argument [13]. 
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The diagrams for operator then give 



3^. 



(1) 



-/ 1 



-/ f 1 



2 
3 



■ 5^x,p - 4J5f,A + 18 J^x,r - 4J^x,y - 5 J^xj + 



2(i82 + i83) 



(3.5) 

^XA - J^x.v 



+a^5^[<^l({Mg)};{M.})^ " I ^^lij ({<};{Mv})^,■ 



P.2]^r,^(7)l._ 



dm 

dJ^XA 



^1 ;•£{<>} 



fl 16 



I ^;-i;({<};{ma})^.v 



(3.6) 



L Dgl({MlJ)};{Mz})J^.:^; 



7e{Mf' } 



(3.7) 



-'^l*B;r J 1 



A' i 16 



(3.8) 



The chiral logarithm functions and c^f'^ are defined in Ref [13], the subscripts on these func- 
tions give the flavor and taste of the relevant meson, S runs over the sea quarks, p sums over taste 
representations P,A,T,V,I (with degeneracies Np), X is the xx meson, and other notation is given 
in Refs. [7, 12]. The contributions to operators G^-G^ have similar forms. 



4. Discussion 

Wrong spin/taste operators arise from the local operator construction and ai^e ^(1) in the lat- 
tice spacing. Their contributions to matrix elements ai^e suppressed to NLO because taste-symmetry 
violation is required. This is different from mixings due to perturbative corrections, which are sup- 
pressed by a^jAn. One-loop in XVY then makes them ^{a^as/An), which is effectively NNLO. 
So non-analytic chiral logarithms do not arise at NLO from the perturbative corrections. A similar 
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statement is true for possible ^ (a) and higher terms in the relation between the naive quarks and 
the spin-taste basis. 

The wrong-spin effects induce mixing in the ZPT of various operators. However, as long as 
all five are analyzed simultaneously, there are no new low-energy constants induced by these 
effects: the jS,, and are all already present in the continuum. 



References [3, 4] focused on the calculation of the quantity t, = (/g, y Bb, ) / {fB,, y Bb^ ) , which 
comes from the matrix element of operator The chiral effects of the wrong spins were not 
known at the time of the HPQCD calculation [3] and were therefore omitted from the analysis and 
error estimate. In the Fermilab/MILC calculation [4], the full ;|^PT expression was available, but 
the matrix elements of operators other than were not calculated, preventing a direct inclusion of 
the wrong-spin effects. However, it was possible to estimate the error of omitting these effects by 
using a small subset of new data to investigate the other matrix elements. The result, § = 1 .268(63) 
included a 3.2% error from this effect, which was the second largest source of error. In the ongoing 
second-generation Fermilab/MILC project [14], matrix elements of all five operators i^,, are being 
computed, which means that the full ;|^PT expression can be used in the analysis, and there will be 
no "wrong-spin error." Of course, a chiral/continuum extrapolation error will remain. 

I thank J. Laiho, R.S. Van de Water, and C. Bouchard for discussions and for help with various 
aspects of the calculation. This work has been partially supported by the Department of Energy, 
under grant number DE-FG02-91ER40628. 
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